Current-current correlations in a model CuO^ system 
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We study a 3-band extended Hubbard model on a CuOz system by the density matrix renormal- 
ization group (DMRG) method. Our system has geometry similar to that of a section of the CuOi 
plane of the copper-oxide superconductors. We have studied the effects of nearest-neighbour repul- 
sions on the current-current correlations of two model CuOz systems. We have calculated several 
types of staggered current correlations as a function of various model parameters such as interac- 
tion strengths and filling. We show that repulsive interactions are clearly capable of increasing the 
magnitude of the current-current correlations in both models studied. The long-distance behavior 
of the current-current correlations is qualitatively influenced by V. We show that repulsive interac- 
tions are capable of enhancing correlations involving at least one copper-oxygen rung over distances 
comparable to the length of the system. In some cases, V decreases the exponent corresponding 
to the power-law decay, leading to a considerable slowing down of the long-distance decay of the 
current-current correlations, compared to that of the Hubbard model. 

PACS numbers: 71.10.Fd, 71.10Hf, 74.20Mn 



I. INTRODUCTION 

The unusual normal-state phenomenology of the high 
T c superconductors constitutes an important theoreti- 
cal challenge since there is not yet a consensus in the 
community about a consistent microscopic theory that 
would reproduce all aspects of the phase diagram of 
the cuprates. Angle-resolved photoemission spectroscopy 
(ARPES) experiments have revealed several properties of 
the so-called pseudogap phase in the underdoped part of 
the phase diagram jjj and it is likely that this phase holds 
the key to a general understanding of the phase diagram. 
Among the scenarios proposed to explain the pseudo- 
gap phase, we mention the 'circulating-current' order- 
parameter, which breaks time-reversal symmetry and 
four- fold rotational symmetry about the copper atoms, 
while preserving the translational symmetry of the lattice 
This theory also gives rise to non- Fermi liquid be- 
haviour. The 'd-density wave' (DDW) order-parameter 
based on orbital currents || has been proposed as a 
hidden order parameter for the phase diagram of the 
cuprates. In Fig. 1 we present a schematic picture 
of these circulating currents in the CuO-2. plane of the 
copper-oxides. 

Phases similar to those with DDW or circulating- 
currents order-parameters have been discussed previously 
in the literature under different guises. The so-called or- 
bital antiferromagnet is a phase with broken translational 
and time-reversal symmetries, described for the excitonic 
insulator H and in the two-dimensional (2D) Hubbard 
model ||. 'Flux phases' were studied in the context 
of the large n limit of the Heisenberg model ||. Stag- 
gered currents were found to be stabilized by attractive 
nearest-neighbour interactions in the 2D Hubbard model 
. On the other hand, t — J clusters were shown 
variationally and numerically to have a staggered flux 



order-parameter. Recently, the doped two-leg t — J lad- 
der, studied by the density matrix renormalization group 
(DMRG) method [[l0| was found to have exponentially 
decaying rung-rung current-current correlations. How- 
ever, the staggered flux phase was found to be stabi- 
lized in the half-filled 2-leg ladder by abelian bosoniza- 
tion methods Jtl| . Further, inclusion of longer-range in- 
teractions in the 2-leg ladder (extended t — U — V — J 
model) appears to stabilize the staggered flux phase, as 
seen from DMRG calculations fll2] . 
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FIG. 1. Schematic picture of the circulating-current phase 
in the copper-oxides. The filled circles represent copper atoms 
and the vertices oxygen atoms of the CuOi plane. The arrows 
on each bond correspond to the direction of the current along 
the bond. 

With this collection of stimulating results on flux 
phases in the Hubbard and t — J models, it would be in- 
teresting to study a system whose geometry more closely 
resembles that of the copper-oxides. This is especially 
relevant in the light of the debate regarding the suitabil- 
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ity of the t — J model to the copper-oxide system and in 
particular the importance of an explicit treatment of the 
oxygen atoms [|l3| . There is now general agreement in the 
community that the CuOi planes play a crucial role in 
the high T c superconductors Q and that the geometry of 
the copper-oxides is relevant to superconductivity. Fur- 
ther, several models with effective interactions have been 
proposed to describe the CuOi plane |l4| , pl||l3| ] . These 
studies have emphasized the role of strong electronic cor- 
relations. Mean-field studies |2| have indicated that it is 
necessary to go beyond on-site interactions and an ex- 
plicit treatment of the nearest-neighbour interactions is 
necessary to stabilize the circulating-current phase. An 
accurate treatment of the strongly correlated full 2D 
CuOi system, though desirable, is not presently forth- 
coming. Therefore, we have opted for a C11O3 system 
that we believe to incorporate the crucial geometric in- 
gredients of the system, including an explicit treatment 
of oxygen sites, which can then be studied by the highly 
accurate DMRG method Jl§ . 




FIG. 2. (a) Single unit cell of the C11O3 system studied, 
the copper site is marked by a filled circle and all vertices 
correspond to oxygen sites. The phase conventions of the 
transfer integrals, corresponding to the symmetry of the d- 
and p-orbitals are indicated (b) Model 1 : 6 unit cells of the 
extended CuOs system with t\\ = (c) Model 2 : similar to 
(b) with t\\ / 0. 

In this paper we study currents in a 3-band extended 
Hubbard Hamiltonian on a Cu(93-type system (Fig. 2), 
which is believed to model the essential aspects of the 
2D system. It is possible to envisage several quasi-lD 
sections of the 2D C11O2 plane. It is necessary to choose 
geometries that are capable of supporting circulating cur- 
rents. The plaquette that forms a single unit cell of the 
CuO-2 plane is shown in Fig. 2a. The system shown in 
Fig. 2b, henceforth to be called Model 1, consists of such 
plaquettes connected at a single point (tn = 0). A strip 
of the CuOi plane is shown in Fig. 2c. In this system, 



henceforth known as Model 2, links parametrized by in 
have been added between adjacent plaquettes to model 
the 2D nature of the system. This is expected to better 
correspond to scenario in the C11O2 planes where the pla- 
quettes are all coupled due to hopping along both axes. 

The essential difference between models 1 and 2 is that 
since the plaquettes of Model 1 are connected at just a 
single point, it is possible to globally reverse the signs of 
the currents in each unit-cell completely independently of 
all the other unit-cells, without any energy cost. There- 
fore, while each unit-cell is a-priori capable of sustain- 
ing currents, long range order is not possible in Model 
1. Thus, it would be necessary to include links between 
these plaquettes (present in the physical system), a sce- 
nario which corresponds to Model 2. In this paper, we 
have examined both these models in detail in order to 
study the effects of interactions locally in a single unit 
cell as well as the long-range effects brought about by 
coupling plaquettes as in Model 2. 

For weak interactions, ladder-type systems can be 
treated using a combination of bosonization and pertur- 
bative renormalization group techniques, as in the case 
of the 2-leg ladder (with disorder) |l7j or N-leg Hubbard 
ladder |Q. In fact, abelian bosonization methods have 
been used to study staggered currents in the 2-leg ladder 
PJ. These powerful analytical tools have proved useful 
in extracting a great deal of information. However, these 
methods cannot fully take into account the exact values 
of the parameters as well as the precise geometry of the 
system, both of which are crucial to the copper-oxides as 
discussed above. While our system appears to be a quasi- 
lD system, it has a highly non-trivial geometry and up 
to 4 th neighbours are connected. Given the large number 
of parameters in the 3-band Hubbard model and the im- 
portance of the delicate balance between them, we have 
chosen to use the highly accurate DMRG method which 
gives information even in the strong-coupling regime in- 
accessible to perturbative approaches. 

We have used a 3-band extended Hubbard model on 
the systems in Fig. 2, described in greater detail below. 
The 3-band Hubbard model for clusters of the CuOi sys- 
tem (up to 6 x 6 unit-cells) has been studied by the con- 
strained path quantum Monte Carlo method (CPMC) 
However, the focus of the CPMC studies was 
pair-binding of holes as well as the pairing and other cor- 
relation functions, but not currents. Thus, a study of 
flux- type phases or currents in this model does not exist, 
to our knowledge. 

II. MODEL AND METHOD 

The 3-band Hamiltonian for our CuOs system in the 
hole picture is given by: 

H = H + H 1 +H 2 (1) 
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Hq — ~tpdd ia (p x ,i+ a /2,cr — Px,i-a/2,a ~ Py,i+a/2,cr 

i , (7 

+ Py,i-a/2,a) ~ tppP Xti+a / 2ja (~Py,i+a/^ + Py,i-a/2,tr) ~ 
{Py,i+a/2,<r ~ Py,i-a/2,<r) 

+ ^2 t \\(pi,i+a/2,aPyJ+a/2,° + Pi,i-a/2,aPy,3~a/2,°) 

+ h.c. + ^2^ P n l Px(y) ^ + e d n da ) 

H 2 = V pd + <,„ ) + Vpp E n L,Ay„ 

i^i ip 

+ Vn N n l n J 

where creates a hole of spin cr on the copper c?-orbital 
at unit-cell i, p x (y),i±a/2,a creates a hole of spin a at the 
oxygen p x (p y ) orbital at sites fi = i± a/2 in the x and y 
directions of unit-cell i. Here (ij) corresponds to a sum 
over nearest-neighbour unit cells i and j and the sum- 
mation over fi corresponds to z ± a/2. Operator n l da is 
the number operator for the copper atom of unit cell i 
and n p ( ^ a are the number operators for the p x (p y ) 
orbitals at sites fi = i ± a/2 of unit cell i. The param- 
eters e d and e p are the site-energies of the copper and 
oxygen sites (A = ed — e P )- The copper-oxygen hopping 
is parametrized by t pd and the oxygen-oxygen hopping 
by t pp while t\\ is the hopping between the p y orbitals 
of neighbouring unit cells. U d and U p are the on-site 
Hubbard repulsions at the copper and oxygen sites, re- 
spectively. The nearest-neighbour copper-oxygen repul- 
sion is parametrized by V pd , the oxygen-oxygen repul- 
sions within a unit cell by V pp and the repulsions between 
p y orbitals in neighbouring unit cells by Vji . Further, the 
phases of the hopping integrals t, corresponding to the 
hybridization of d- and p-orbitals are represented in Fig. 
2. 

We study this model by the DMRG method. The com- 
plicated geometry (up to 4th neighbours are connected) 
of this system required the implementation of a 7-block 
scheme (Fig. 3). This choice of RG scheme ensures that 
the relatively complicated geometry, believed to be cru- 
cial in the cuprates, is preserved at each step of the calcu- 
lation. In Fig. 3a, we see the breakup of the two unit-cell 
system which is the starting point of our calculation. The 
system size increases by one unit cell (4 sites) at each it- 
eration. The blocks marked T' and '7' in the figure are 
renormalizcd at each iteration and five new blocks (un- 
renormalized) of one site each are added to the middle 
of the system. The dimension of the Hilbert space of 
each of the added sites is 4, corresponding to empty, one 
hole (up and down spins) and 2 holes. In Fig. 3b, we 



represent the blocks and superblocks of the three unit- 
cell system. The following iterations can be obtained by 
replacing the superblocks T' and '7' by the appropriate 
system. The DMRG method is a highly efficient way 
to study 1- and quasi 1-dimensional strongly correlated 
fermionic systems. We are also not limited to a choice of 
fillings that correspond to closed-shells, a habitual con- 
straint in quantum Monte Carlo simulations. 



(a) 




FIG. 3. The breakup of the extended system into blocks 
and superblocks used in our 7-block DMRG scheme for (a) 
even iterations and (b) odd iterations. Copper atoms are 
marked by filled circles and oxygen atoms by open circles. The 
physical bonds are marked in dashed lines. The superblocks 
of each iteration are marked in dot-dashed lines and also carry 
their number (1 and 7). The five central blocks consist of one 
site each. Superblock sizes increase by two sites in each iter- 
ation and the superblock of the following iteration is marked 
by solid lines. 

The parameters varied in our study are mainly the 
interaction strengths (V's), A and the filling. We have 
studied this model with t pp = 0.6, U d — 10.0 and U p = 
5.0 in units where t pd = 1.0. Further, we consider two 
different geometries corresponding to values of i|| = 
and t\\ = 0.6. The former corresponds to the case of 
plaquettes that are coupled only at a point (Model 1) 
and the latter to the case of fully connected plaquettes 
(Model 2). The choice of parameters for 3-band models 
of the copper-oxides has been extensively discussed in the 
literature . Our choice corresponds to quite a realistic 
range of parameters for the copper oxides []ig| — |21f| . In fact, 
one of the inherent difficulties associated with any study 
of the copper-oxides is the large number of parameters, 
giving rise to a huge phase space. We have tried to vary 
the different V's in order to locate the interesting part 
of the phase diagram. One of the difficulties encountered 
in our calculations has been the inaccessibility of certain 
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ranges of V due to convergence problems, as discussed in 
the end of Section IIIB. 

We have studied up to 25 CuO^ unit cells, correspond- 
ing to 101 sites, at different dopings. Half-filling corre- 
sponds to one hole per copper site in the hole picture, 
where the system is known to be an antifcrromagnetic 
insulator for a wide range of parameter values. We have 
studied systems with different numbers of holes, corre- 
sponding to doping of the physical system. We add two 
holes at every other iteration. Thus, the filling is not 
the same at every step of the calculation. The cases we 
discuss are systems with 20 - 25 unit cells. We have 
studied systems of N unit cells with holes given by, 
rih = N + 3 + mod(N + 1,2). The filling x is given by 
rih/N — 1. Thus, x = corresponds to the half- filled 
case. We have studied fillings x = 0.12, 0.13 and 0.19 
(N = 25, 23, 21) for odd numbers of unit cells and x 
= 0.17, 0.18 and 0.25 (N = 24,22,20) for even num- 
bers of unit cells. Thus, we are able to study a range of 
doping values. In addition, we have studied the case of 
n h = N + 1 + mod(N + 1, 2), for the model with t\\ = 0. 

We compared the ground state energies for several it- 
erations with the energies obtained from exact diagonal- 
ization calculations to verify our RG scheme. We var- 
ied the number of states m on the superblocks ( m = 
125, 150 and 200). We have kept up to 200 states on the 
superblocks. This correspond to a Fock space dimension 
H F = m 2 x 4 5 = 40,960,000. This provides sufficient ac- 
curacy, estimated in several ways. The discarded weight 
of the density matrix states is of order 10 -6 with 200 
states. In certain cases, we have also calculated the en- 
ergies and properties of the non-interacting system (free- 
particle model with all U's and V's set to zero), where 
exact results are available. This provides a strong check 
on our data, since convergence is known to be most diffi- 
cult in this case for the DMRG algorithm. Another prop- 
erty that we have used to check the quality of our data 
is the magnetic moment | (£?*)!, at site i, which ought to 
be zero in principle. Thus, only data which correspond 
to relatively small magnitudes of \{S?)\ have been used. 

III. RESULTS AND DISCUSSION 

The physical property we have studied is the current- 
current (CC) correlation function for links k and I. The 
current along a bond connecting sites [i and v is de- 
fined as = it^J^^la^o- ~ 4a c A»o-)' Tne mean - 

field analysis of the 3-band Hubbard model indicated that 
the nearest-neighbour repulsion V is essential to stabilize 
the circulating-current phase || . The circulating-current 
phase shown in Fig. 1 would then consist of currents rep- 
resented by the arrows in the figure. Since we work in 
the real representation of the Hamiltonian, we use the 
current-current correlation functions as a measure of this 



order. 

We compared our data for the non-interacting model 
with the exact data available in this case to assure our- 
selves of convergence. The current for bond I in unit-cell 
i is denoted by J % {1), according to the conventions for 
bond numbers shown in Fig. 4. We refer to the vertical 
copper-oxygen bonds as rungs and the horizontal copper- 
oxygen bonds as legs in analogy to the terminology used 
for ladders. 



(a) 




i=(N-lV2 i=0 i=0' 



FIG. 4. Conventions for current- current correlation func- 
tions studied for systems with (a) odd and (b) even numbers 
of unit cells. Unit-cell numbers i are marked beneath and the 
numbers in the figures correspond to the bond numbers. 

We have calculated the equal-time, zero-frequency 
current-current (CC) correlation function, which is the 
ground state expectation value (J°(k)J l (l)). These cor- 
relations are between currents in the central unit cell, 
denoted by J°(k) and currents on the other unit-cells on 
one half of the system J z (l) (i = 1, . . . (N - l)/2), for 
systems with odd numbers of unit-cells. Slightly differ- 
ent correlations are calculated at even and odd iterations 
due to the DMRG scheme used and these are pictured 
in Fig. 4. For systems with even numbers of unit-cells, 
we introduce the notation (J (k),P(l)), and ( J ' (fc) J*(Z)) 
with i — 1, . . . (N — l)/2, as shown in Fig. 4. However, 
we continue to refer to both cases as (J°(k)J l (l)) in the 
text, with the different notations for even and odd unit- 
cells being implicit. By the bilateral symmetry about the 
O — Cu—O axis, it is not necessary to take into account all 
the possible cases. We therefore restrict ourselves to the 
case of (J°{k)J l {l)), k,l — 1, . . .3. We analyse cases of 
systems with 20-25 unit-cells, since these are sufficiently 
large to study the decay of the current-current correla- 
tions. In particular, we analyze the behaviour of the dif- 
ferent types of diagonal ((J°(k)J l (k)) and off-diagonal 
{J°(k)J l (l)), (k =fi I) current-current correlations and try 
to classify them. 
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A. Analysis of Model 1 

We first present our results for the CC correlations on 
the CuOz system with Model 1 [t\\ = 0) which consists 
of plaquettes coupled at a single point. We have per- 
formed calculations with values of t p d and t pp equal to 
1.0 and 0.6. We have varied A and V p d, while V pp was 
set to 0. We have further considered two sets of fill- 
ings corresponding to rih = N + 3 + mod(N + 1,2) and 
rih = N+ 1 + mod(N+ 1, 2). The second series then gives 
us fillings x = 0.04, 0.043 and 0.048 (N = 25, 23, 21) and 
x = 0.083, 0.09 and 0.1 (N = 24,22,20). Thus, the 
second series rih = N + 1 + mod(N + 1,2) has the ad- 
vantage of providing fillings that are practically constant 
for different system sizes. This helps to verify that our 
conclusions hold over different system sizes. 
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FIG. 5. Current-current correlations |{J°(fc) J l (Z))j for a 
25 unit-cell system with 28 holes, for Model 1, in a log- log 
representation. All systems have the parameters t pp — 0.6, 
tu = 0, A = 1 in units with t v a — 1. Circles represent the 
case with the above parameters and all V's = 0, squares are 
for Vpd = 1 and diamonds are for V p d = 2, all with Ud = 10 
and Up — 5, obtained from DMRG calculations with m = 200 
states. Triangles represent the exact values for the correlation 
for the non- interacting model (U p = Ud — V p d = V pp = 0), 
and inverted triangles the values obtained for the same pa- 
rameters from m = 200 state DMRG calculations, presented 
for comparison. Dashed lines are a guide to the eye while 
full lines correspond to the power-law fits of the data, (a) 
represents the [(J°(2)J i (l))| correlation with a* NI = -2.94, 
ckni = —3.12, ao = —3.09, a\ = —3.03, a.i = —3.09 
and (b) the |(J°(2)J i (3))| correlation with a* NI = -2.92, 
oni = —3.11, ao = —3.16, ai = —3.21, «2 = —3.37. See 
text for definitions of the as. 

We have also carried out numerical fits on our data. 
For our results on Model 1, we found our data well- fitted 



by power-laws. We define a^i as the slope corresponding 
to the non-interacting case obtained from DMRG and 
a* NI as the slope of the exact data. Likewise at,k = 
. . . 2 correspond to the cases with V p d = k. These slopes 
give us the exponents of the power-law. 

We have studied the effect of increasing the copper- 
oxygen repulsion V p d on the current-current correlations. 
In Fig. 4a, we present results for the dependence of the 
CC correlation function \ {J°(2) J l (l))\ with distance i for 
a 25 unit-cell system with 28 holes. From Fig. 4, it can 
be seen that this is the correlation of a rung of the central 
unit-cell with an oxygen-oxygen bond (2 = 1) of the unit- 
cells away from the centre. We first consider this corre- 
lation for the non-interacting case (all t/'s and Vs set to 
0), where exact results are available. We have compared 
the results obtained for |(J°(2)J i (l))| from a DMRG cal- 
culation with to = 200 states with the exact results. The 
convergence of the data is excellent, as seen in Fig. 5a. 
The exponent of the DMRG data a^i — —3.12, com- 
pares favorably with the exponent obtained for the exact 
data, a* NI — —2.94. This is a strong check on our results, 
since the convergence of the DMRG method is known to 
be harder for the free-particle model than for the equiv- 
alent system with strong correlations. 

We now note from Fig. 5a that the curve which rep- 
resents current-current correlations for the purely local 
Hubbard model with Ud = 10, U p = 5 and Vs = 0, 
lies above the non-interacting curve. We then turn on 
nearest-neighbour interactions. Increasing the copper- 
oxygen repulsion V p d = 1 has the effect of increasing 
the magnitudes of the current-current correlations. Fur- 
ther increasing V p d to 2 increases the correlations, but 
by a smaller amount. The decay of the correlations fol- 
lows power-law behaviour, as seen from the log-log plots. 
Thus, we observe that repulsive interactions are capable 
of increasing the magnitude of the currents. It can be 
seen from the curves that the decay at large-distances is 
hardly affected by the Vs, the curves are nearly parallel. 
This is borne out by the slopes (—3.09, —3.03 and —3.09), 
which are practically unchanged. This is because the pla- 
quettes are coupled at a single point, which prevents the 
formation of a state with long-range-order. 

In Fig. 5b, we present our results for the correlation 
\(J (2)J'(3))|. This is a correlation between a copper- 
oxygen rung of the central unit-cell, with a copper-oxygen 
leg of the side unit-cells. Increasing the copper-oxygen re- 
pulsion enhances the magnitude of the correlations com- 
pared to the Hubbard model, while their decay is power- 
law, again with similar exponents, as is seen from the 
curves which are nearly parallel and from the power-law 
fits. 

In Fig. 6a we present the behaviour of the correlation 
|(J°(2) J'(2))|, which is the copper-oxygen rung-rung cor- 
relation function. The data presented are for the case of 
25 unit-cells with 26 holes, since the absolute magnitude 
of the correlations is much larger than in the case of 28 
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holes. The qualitative behaviour remains the same. We 
observe that the absolute values of the V p d = 2 data are 
enhanced by almost an order of magnitude compared to 
the Vpd = curve. We have fitted our data with two se- 
ries of power-laws for even and odd sites in order to take 
into account the alternating behavior of the curves. The 
slight differences between ao and aj can be attributed to 
the rather large alternance of the V p d — data. This al- 
ternance is reminescent of the generic behavior of correla- 
tion functions in quasi-lD systems, due to the harmonics 
of hp. We notice that the magnitude of the oscillations 
decreases when V is increased. In quasi- ID systems, sim- 
ilar behavior of spin-spin and charge-charge correlations 
is interpreted as an increase in the value of the charge- 
stiffness (K p ) HJ. 
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FIG. 6. Current-current correlations for a 25 unit-cell sys- 
tem for Model 1. Note the log-log presentation. All systems 
have the parameters (t pp — 0.6, t|| = 0, A = 1, Ud = 10, 
Up — 5). Circles represent the case with the above parame- 
ters and all V's = 0, squares are for V p d = 1 and diamonds 
are for V p d = 2. Dashed lines are a guide to the eye, while 
full lines represent the power-law fits, (a) |(J°(2)J ! (2))| for 
a 25 unit-cell system with 26 holes with (ao = —3.34 (even) 
and -3.94 (odd)); (e*i = -4.00 (even) and -4.18 (odd)), and 
(b) |(J°(3)J i (3))| for a 25 unit-cell system with 28 holes with 
(q = -2.40, ai = -2.69, a 2 = -3.00). 

The behaviour of the correlations |(J°(3) J l (3))| can be 
seen in Fig. 6b. This correlation is a copper-oxygen leg- 
leg correlation function, as seen from Fig. 4. Here, the 
behaviour observed is different from the previous cases 
studied, since the correlation is enhanced very slightly 
by increasing repulsions for first neighbours. However, 
the situation is reversed at very short distances and in- 
creasing repulsions causes the correlations to decay even 
faster, as can be noted from the crossing of the curves 



in Fig. 6b. Their decay is power-law with the exponents 
decreasing, as can be seen from the figure. 

We have also analysed the behaviour of correlation 
functions of the type | (J°(1)J 4 (1)) |. These are not 
presented for reasons of space. These correlations are 
oxygen-oxygen correlations and decay as power-laws, 
with the magnitudes falling as V p d is increased. This 
however, could be due to the value of V pp = 0, chosen 
here. 

We now briefly discuss the case of x — 0.18 (22 unit- 
cells with 26 holes), which fits into the above analysis 
with an interesting reversal. In this case, the rung-rung 
correlations (indeed, correlations involving one rung) are 
initially decreased when V p d is turned on. Increasing V p d 
to 2 brings about a reversal of this behaviour, with the 
magnitude of the correlations increasing compared to the 
V p d = case. In Fig. 7 we present one example of this 
behaviour with the correlation \ (J° (2)J I (3)}|, which is a 
rung-leg correlation. The reversal of the trend is clearly 
seen from the V p d = 1 curve (lies below the V p d = curve) 
and the V p d — 2 data which is enhanced with respect to 
the V p d = curve. This indicates that there could be a 
threshold value of nearest-neighbour interactions V for 
enhancing the magnitude of current-current correlations. 
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FIG. 7. Current-current correlations |(J°'(2)J ! (3))| for a 
22 unit-cell system with 24 holes. Note the log-log represen- 
tation. All systems have the parameters (t pp — 0.6, iy = 0, 
A = 1, Ud = 10, Up = 5) Circles represent the case with the 
above parameters and all V"s = 0, squares are for Vpd = 1, 
diamonds for V p d = 2. Lines represent power-law fits with 
exponents ao = —2.41, ct\ — —2.69 and a 2 = —3.00. 

Thus, the above analyses show that the copper-oxygen 
rung-rung correlation (or even correlations involving one 
rung) are enhanced in magnitude by repulsive interac- 
tions, compared to the purely local Hubbard model, 
probably requiring a minimum threshold value of V. 
The copper-oxygen leg-leg correlations show an initial in- 
crease followed by a decay in magnitude and correlation 
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length, as V p d is increased. The oxygen-oxygen correla- 
tion is diminished in magnitude and correlation length, as 
Vpd is increased. We have analysed the cases of different 
fillings described above and have found similar behaviour 
in all these cases. 

B. Analysis of Model 2 

We now turn to an analysis of Model 2 with tu = 0.6, 
The behavior of this model is expected to be quite differ- 
ent from that of the tu = case. It is expected that the 
long-distance behavior of the CC correlations is much 
more sensitive to the nearest-neighbour repulsions V 
Here, we first study the purely local Hubbard model on 
our lattice (Ud = 10, U p = 5 and Vs = 0). We note that 
this corresponds to the system in Fig. 2c, with bonds 
connecting different plaquettes (t\\ = 0.6). However, the 
V on this bond is initially set to zero, to give us a model of 
connected plaquettes with purely local interactions. We 
the begin to turn on Vs, i.e we turn on copper-oxygen 
and oxygen-oxygen nearest-neighbour interactions, along 
with the interaction Vj| ( V p d — V pp = Vj| = 1.0). There- 
after, we further increase V\\ to study the effect of longer- 
range interactions on the current-current correlations. 

We shall first consider the case of 25 unit cells at a 
filling of 0.12 in some detail before analysing our other 
data. In Fig. 8a - 8c we present diagonal correlations 
of the type \{J°{k) J\k))\, for k = 1, . . . 3, on a log- 
log scale. These correspond respectively to the oxygen- 
oxygen diagonal correlation, the rung-rung correlation 
and the leg-leg correlation of bonds in the central unit 
cell with bonds in the side unit cells. In Fig. 8a, we 
present the decay of the oxygen-oxygen diagonal correla- 
tion (|(J°(1) J l (l))|) with increasing distance for increas- 
ing values of the nearest-neighbor repulsion V These 
are compared to the data with all Vs = 0, i.e. the 
purely local Hubbard model. From the data in Fig. 8a, 
we note the effect of repulsions on the correlations be- 
tween currents on oxygen-oxygen bonds. Turning on Vs, 
(Vpd — Vpp — Vj| =1) causes the correlations to decay 
in magnitude compared to the Hubbard model. Further 
increasing Vji to 3, causes the correlations to decay fur- 
ther. Data for V p d = V pp = 1, Vj| = 2, fit smoothly into 
this pattern and have not been shown, for clarity. The 
correlations show exponential decay as can be seen from 
the plots as well as the inset. We have fitted the Hub- 
bard model curve and the Vji = 3 data (scaled by a factor 
of 10) by expressions of the form y = Aexp(px), shown 
in the inset. The clear straight lines seen in the semilog 
plot are indicative of the quality of the fit. We obtain 
A = -6.0065 and f3 = -0.6296 for the Hubbard model 
(Vs = 0) and A = -3.3280 and (3 = -0.7753 for the 
Vj| = 3 data. 
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FIG. 8. Diagonal current current correlations for a 25 
unit-cell system with 28 holes (|(J°(fc) J l (k))\) in a log- log rep- 
resentation. All systems have the parameters (t pp = t\\ = 0.6, 
A = 1, Ud = 10, Up = 5) in units with t p d = 1. Cir- 
cles represent the case with the above parameters and all 
Vs = 0, squares are for V p d = V pp = Vy = 1, triangles 
for V p d = V pp = 1, Vj| = 3. Triangles are in some cases 
filled, for clarity. Lines in the main figures are meant to be 
a guide to the eye. Lines in the insets correspond to fits 
which are discussed in greater detail in the text. From top 
to bottom, (a) diagonal oxygen-oxygen correlation function 
|(J°(1)J 1 (1))|. Inset (semilog representation) contains V\\ = 
3 data multiplied by a factor 10, for visibility, (b) rung-rung 
correlation |(J°(2) J*(2))|. The inset (log-log plot) contains 
Vji = 3 data multiplied by a factor 100. (c) leg-leg correlation 
|(J°(3)J < (3))|. The inset (log-log plot) contains Vy = 3 data 
multiplied by a factor 10. 

In Fig. 8b we present correlations of type 
|(J°(2) J l (2))|, which correspond to the copper-oxygen 
rung-rung correlation, for different values of interaction 
strengths. Here, we notice a sizeable enhancement of 
the current-current correlations (by an order of magni- 
tude) over distances of a few lattice spacings as inter- 
actions (Vs) are turned on, compared to the Hubbard 
model. Thereafter, interactions cause the correlations to 
fall, with respect to the Vs = case. However, further 
increasing Vji causes the correlations to rise i.e. there is 
a reversal of the trend. It is reassuring that the DMRG 
data in these cases are well converged, as seen from the 
small discarded weight and the small magnetic moments 
at the central sites, \Sf\ < 4xl0~ 4 . These data have been 
fitted by power-laws. The fits are good as seen from the 
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straight-lines in the log-log plots, both in the main figure 
and in the inset (thick lines correspond to the fits). The 
slope of the straight line gives the exponent of the power- 
law and we obtain a Q = -4.4070 and a 3 = -4.6830. We 
also note that the cross-rung-rung correlation, of type 
( J°(2) J l (4)) | are similar qualitatively and quantitatively 
to the rung-rung correlations and are not presented here. 

In Fig. 8c we present correlations of type 
(J°(3)J l (3))|, which correspond to the copper-oxygen 
leg-leg correlation, for different values of interaction 
strengths. Here, we notice a slight enhancement of the 
current-current correlations over distances of a few lat- 
tice spacings as interactions (V's) are turned on, com- 
pared to the Hubbard model. Thereafter, interactions 
cause the correlations to fall, with respect to the V's = 
case, upon increasing Vs. The long-distance behaviour 
in this case is interesting. The V = data (Hubbard 
model) are clearly exponential as seen from the main fig- 
ure and the inset (A = -3.9204 and j3 = -0.6254). On 
the other hand the V\\ = 3 data (triangles) show much 
more complex behaviour corresponding to power-law be- 
havior at short distances (up to 6 sites) followed by an 
exponential tail at long distances. The power-law fit is 
indicated in the inset by the thin line and the exponential 
tail by the thick line. We obtain the power-law exponent 
a 3 = —2.6299 from the slope of the power-law part of the 
fit and A = -0.0143 and = -0.8947 for the exponential 
part (thick line). 

We now turn to an analysis of the off-diagonal cor- 
relations. Given the large number of correlations, we 
attempt to set out the general trends by treating one off- 
diagonal correlation of each type. In Fig. 9a - Fig. 9c, 
we present off-diagonal CC correlations, |(J°(1) J*(2))|, 
\(J°(3)J i (2))\ and ( J°(3) J l (l)) |. The first of these cor- 
responds to the correlation of an oxygen-oxygen bond 
in the central unit-cell with copper-oxygen rungs of the 
side unit cells. We refer to this henceforth as an "oxygen- 
rung" correlation. As seen from the log-log plot of Fig. 
9a, the correlation (1,2) exhibits power-law decay both 
for the Hubbard model and upon increasing Vs. There 
is a substantial increase of this correlation at short dis- 
tances with increasing V. This trend persists over long 
distances for the case of Vj| =3 as seen from the fig- 
ure. The power-law fits to these data are presented in 
the inset, with the V = 3 data scaled by a factor 10. We 
have fitted the Hubbard model data in two series corre- 
sponding to the odd and even points to take care of the 
alternance present in our data and obtain aos of —3.9408 
(odd) and —3.3321 (even). The Vfi = 3 data were fitted 
by a power-law with exponent a 3 = —4.0923. We note 
that the alternance of the V\\ — 3 is strongly diminished 
with respect to the V = data. 

Similar behavior is observed in the leg-rung correlation 
function |(J°(3) J 4 (2))|, presented in Fig. 9b. The V\\ = 3 
data are clearly enhanced with respect to the Hubbard 
model data almost over all the length of the curve. We 



followed a procedure similar to the previous case, fitting 
the V = data by power-laws in two series to obtain 
a s of -3.0559 (odd) and -4.0292 (even). The power- 
law fits to the Vji = 3 data gave us an exponent 0:3 of 
—4.1966. Correlations quench the even-odd alternation 
of the data. 
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FIG. 9. Off-diagonal current current correlations 

|(J°(fc)J 4 (/))| (k I) for a 25 unit-cell system with 28 holes 
in a log-log representation. All systems have the parameters 
(t pp = t\\ — 0.6, A = 1, Ud = 10, Up — 5) in units with 
t p d — 1. Circles represent the case with the above parame- 
ters and all V's = 0, squares are for Vj| = V pp = Vy = 1, 
triangles for V p d = V pp = l,Vj| = 3. Symbols are in some 
cases filled, for clarity. Lines in the main figures are meant 
to be a guide to the eye. Lines in the insets correspond to 
fits which are discussed in greater detail in the text. From 
top to bottom, (a) "oxygen-rung" correlation |(J°(1)J'(2))|. 
Inset (log-log) contains Vy = 3 data multiplied by a factor 
10, for visibility, (b) leg-rung correlation | ( J°(3) J*(2)> |. The 
inset (log-log plot) contains Vji = 3 data multiplied by a fac- 
tor 10. (c) "leg-oxygen" correlation |(J°(3) J'(l))|. The inset 
(semilog plot) contains Vy = 3 data multiplied by a factor 10. 

In Fig. 9c, we present the correlation |(J°(3) J l (l))|, 
which corresponds to the correlation of the central 
oxygen-oxygen bond with the copper-oxygen leg of the 
side unit cell. We refer to this correlation henceforth as 
the " oxygen-leg" correlation. We observe from the figure 
and from the fits in the inset that these correlations have 
exponential behavior for the Hubbard model as well as 
for the case with Vji = 3. From the fits seen in the inset 
(presented in semilog form) we obtain A = —4.61476 and 
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(3 = -0.6600 for the V — case and A = -1.849094 and 
(3 = -0.8081 for the V\\ = 3 case. 
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FIG. 10. Current-current correlations |( J ' (2) .P (fc)> | for 
a 22 unit-cell system with 26 holes. All systems have the pa- 
rameters t pp = £|| = 0.6, A = 1, Ud — 10, Up — 5 in units with 
t p d = 1. Circles represent the case with the above parameters 
and all V's = 0, squares are for V p d = V pp = Vj| = 1 and 
triangles for V pt j = V pp = 1, Vj| = 3. Triangles are in some 
cases filled, for clarity. Lines in the main figures are meant as 
a guide to the eye and lines in the insets are the power-law 
fits to the data. In both insets, Vy = 1 data (squares) is 
multiplied by a factor of 5 and Vy = 3 data (triangles) by 
a factor of 25 for clarity. From top to bottom, (a) rung-leg 
correlation, |(J°'(2)J < (3))|. From the fits, a = -4.4157, 
ai = —4.3648 and a.i = —3.5904. (b) another rung-leg 
correlation (J ' (2)^(5)) |. From the fits, a = -4.9475, 
ai = -4.7718 and a 2 = -4.3446. 

In addition to the data presented in Figs. 8-9, for 
the 25 unit-cell system with filling x = 0.12, we dis- 
cuss some results for the case of 22 unit-cells with x = 
0.18. In Fig. 10, we present the evolution of two types 
of rung-leg correlations that are similar: |(J° (2)J 8 (3))| 
and |(J°'(2)J 4 (5))|, henceforth referred to as the (2', 3) 
and (2', 5) correlations respectively. In Fig. 10a, we see 
that the V]| = 1 data of the (2', 3) correlation is already 
enhanced with respect to the Hubbard model data over 
intermediate distances. The Vj| = 3 data is clearly en- 
hanced with respect to the Hubbard model data. We 
have fitted these data by means of power-laws, as seen 
from the log-log scale of the curve. The results of the fits 
are seen from the inset. We see that 0:3, the exponent of 
the V|i = 2 curve is visibly reduced to —3.59 compared 



to a® — —4.94 for the Hubbard model and the V\\ = 3 
data falls off much more slowly than the Hubbard model 
data. As seen from the inset, the number of points we 
were able to include in the fit decreases with increasing V , 
as compared to the Hubbard model. It is accepted that 
convergence of the DMRG algorithm is best in the pres- 
ence of strong, local interactions. Our Hubbard model 
case falls into this category with Ud = 10 and U p = 5. 
Increasing V clearly renders convergence much more dif- 
ficult, especially when Vj| — 3 which is most probably the 
reason for this behavior. 

In Fig. 10b, we present in a manner analogous to Fig. 
10a, data for the (2', 5) correlation functions. As in the 
previous case, we note that the Vfj = 3 data is stabilized 
compared to the V = data and the Vj| = 1 data, with 
the exponents ao = —4.9475, u\ = —4.7718 and a$ — 
-4.3446. Thus, we note that the decay of the (2', 5) 
correlation is slower for Vfi = 3 than for the Hubbard 
model. We have dropped 3 points from the Vji = 3 curve 
and 1 point from the Hubbard model curve to obtain 
these results. 

Thus, the analysis of the diagonal and off-diagonal 
CC correlation functions of Model 2 yields the follow- 
ing trends : the nearest-neighbor repulsive interaction 

V is clearly capable of increasing the magnitude of the 
current-current correlations locally and of strongly mod- 
ifying the long-distance behavior. The correlations fall 
into different categories. 

The rung-rung correlation function is enhanced by in- 
creasing Vs and shows power-law decay. In certain cases, 
the exponent of the power-law decreases upon increasing 

V and the correlations fall off slowly as compared to the 
Hubbard model. In fact, this extends to any correlation 
that includes at least one rung, diagonal or off-diagonal. 

The leg-leg correlations show exponential decay in the 
Hubbard model. However, increasing the values of Vs 
and Vj| in particular changes this behaviour - we now 
have a power-law behavior at short distances followed by 
an exponential tail. Thus, we see that Vs are capable 
of locally enhancing the CC correlations and by connect- 
ing the plaqucttes, the long-distance behaviour is greatly 
modified. We believe that the difference between the 
rung-rung and leg-leg correlations is due to our choice of 
model and that this difference should disappear in the 
2D system. 

The oxygen-oxygen correlation function shows expo- 
nential decay with very slight enhancement of the abso- 
lute values for Vj| = 3. This is also true of the "oxygen- 
leg" type correlation functions. This could however be 
due to our choice of parameter values. Clearly, it would 
be interesting to study other parameter ranges for the 
Vs, especially for V p d and V pp . This point is discussed 
further down in this section. 

We have studied various cases involving other parame- 
ter values. Notably, we tried to increase V p d and V pp (V p d 
= 2.0, 3.0, V P p — 2.0, in units of t p d), but encountered 
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convergence problems at the level of the Lanczos diag- 
onalization algorithm. It is quite likely that increased 
degeneracy of the ground state slows down convergence, 
as is known to occur in other cases such as the J\ — J2 
spin chain, at parameter values close to the critical point. 
This is most probably due to the increased derealization 
of the ground state, brought about by larger Vs. To 
overcome this problem, it would be necessary to increase 
the number of states m that we keep on the superblocks. 
However, at our present value of m = 200, we are already 
at the limit of our computational resources. The charge 
densities in our systems show that the charges tend to be 
concentrated on the copper atoms, which possibly com- 
petes with the stabilization of currents. Thus, increasing 
V p d and Vp P is likely to bring about a rehomogenization 
of charge and this could be favourable to the current- 
current correlations. 

At this juncture, it is interesting to compare the mod- 
els with tu = and 0.6. We have observed enhancement 
in the magnitudes of the currents in both models. Thus, 
repulsive nearest-neighbour interactions are clearly capa- 
ble of increasing the magnitudes of currents. The long- 
distance behaviour of correlations is different in the two 
models. Model 1 shows power-law decay of correlations 
and the interaction V has a very small effect on the ex- 
ponents. In Model 2, any correlation involving a rung 
shows power-law behaviour. We have very interesting 
indications that repulsive nearest-neighbour interactions 
can enhance currents over distances comparable to the 
length scale of the system. In certain cases, the exponent 
is considerably reduced and the current-current correla- 
tions stabilized compared to the Hubbard model. 



IV. CONCLUSIONS 

In conclusion, we have studied a 3-band extended Hub- 
bard Hamiltonian on a model C11O3 system, with geom- 
etry similar to that of a section of the CuOi planes of 
the copper-oxides. We have studied two different cases, 
with the parameter t\\ — and 0.6, which correspond re- 
spectively to Model 1 with plaquettes coupled at a single 
point and Model 2, with plaquettes coupled all along the 
length of the system. Our calculations show that nearest- 
neighbour repulsive interactions are clearly capable of en- 
hancing current-current correlations in magnitude and of 
qualitatively changing the long-distance behavior of the 
system. Model 1 shows power-law decay of correlations 
and the interaction V p d has a very small effect on the ex- 
ponents. Increasing the copper-oxygen repulsion V pc i can 
increase the magnitude of the rung-rung current-current 
correlation (or correlations involving a single rung) by 
up to an order of magnitude. In Model 2, the effect of 
introducing nearest-neighbour interactions V in the sys- 
tem is to produce a sizeable increase in the rung-rung 
correlations over distances of a few lattice spacings. Vji 



can enhance correlations involving at least one copper- 
oxygen rung over distances comparable to the length of 
the system. In certain cases, interactions bring about 
a sizeble decrease in the exponent, causing the correla- 
tions to decay slowly with respect to the Hubbard model. 
We have indications that threshold values of the nearest- 
neighbour repulsions are necessary to observe stabiliza- 
tion at long distances. Tuning the copper-oxygen and 
the oxygen-oxygen repulsions would redistribute charge 
more uniformly in the system which is likely to stabilize 
the current-current correlation functions in this system. 
We believe that the enhancement and stabilization of cer- 
tain correlation functions seen in our model constitutes 
a trend that will hold for the 2D system. 
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